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SUMMARY 


Scintillating  communication  signals  generated  by  two  models  are 
compared  with  Wideband  satellite  data  in  a simulated  demodulator  for  phase- 
shift  keyed  (PSK)  signals.  To  effect  this  comparison,  a digital  simulation 
of  a PSK  modem  designed  for  DSf.S  II  is  exercised  for  three  different  input 
signals.  One  input  signal,  which  forms  the  basis  for  the  comparison,  con- 
sists of  detrended  data  from  the  DNA  Wideband  satellite  experiment.  A second 
input  signal  is  obtained  from  the  statistical  two-component  log-normal  plus 
joint-gaussian  model  proposed  by  SRI  International  (SRI I ) . 1 The  initial 
statistics  required  for  this  model  are  obtained  from  analysis  of  the  DNA 
Wideband  satellite  data.  The  final  input  signal  is  obtained  from  a numerical 
multiple  phase-screen  (MPS)  simulation  which  was  generated  to  "match"  the 
initial  Wideband  satellite  experimental  data.  The  purpose  of  the  comparisons 
is  to  test  the  capabilities  of  models  now  in  use  to  represent  the  effects  of 
transionospheric  propagation  in  communication  system  simulations. 

Results  indicate  that  both  the  statistical  model  and  the  MPS  model 
yield  receiver  performance  similar  to  that  obtained  with  the  DNA  Wideband 
satellite  data.  Thus,  at  least  for  this  class  of  receiver  and  for  PSK  sig- 
nals, both  methods  of  generating  simulated  signals  are  useful.  The  MPS 
simulation  has  the  potential  advantage  of  yielding  information  about  iono- 
spheric structure  including  rms  electron-density  fluctuations  and  electron- 
density  power  spectral  density.  However,  the  MPS  simulation  requires 
information  on  the  link  geometry  as  well  as  the  ionospheric  structure  as 
input.  The  major  advantage  of  statistical  models,  such  as  the  combined 
log-normal  and  joint-gaussian  model,  is  the  simplicity  of  generating 
signals  with  such  statistics.  However,  such  models  require  information  on 
a number  of  statistical  parameters. 


maasm  * 

ms  i'.'./c  &:$••-«  |jf 

DOC  E.  • Sect  re?.  Q 

wWfflO'JKiEE  a 

JWTfflCATlOH 


BT 


lCTI8tfltOP»lUBW  m 


CONTENTS 


PAGE 

SUMMARY  1 
ILLUSTRATIONS  3 
SECTION  1 -INTRODUCTION  5 
SECTION  2— PARAMETER  ESTIMATION  7 
SECTION  3—2-COMPONENT  MODEL  12 
SECTION  4— MULTIPLE  PHASE  SCREEN  SIMULATION  25 
SECTION  5— RECEIVER  SIMULATION  33 
SECTION  6— SIMULATION  RESULTS  37 
SECTION  7— CONCLUSIONS  40 
REFERENCES  41 


2 


ILLUSTRATIONS 


FIGURE  PAGE 

1 Severely  scintillating  UHF  signal —(detrended  intensity 
and  phase  from  Wideband  satellite  pass  at  Ancon,  Peru, 

on  16  December  1976).  9 

2 Moderately  scintillating  UHF  signal —(detrended  intensity 
and  phase  from  Wideband  satellite  pass  at  Ancon,  Peru, 

on  16  December  1976).  10 

3 L-band  signal  during  severely  scintillating  UHF  condi- 

tions—(detrended  intensity  and  phase  from  Wideband 
satellite  pass  at  Ancon,  Peru,  on  16  December  1976).  11 

4 Severe  UHF  scinti 1 lation— scatter  component.  14 

5 Severe  UHF  scintil  lation— focus  component.  15 

6 Moderate  UHF  scintillation— scatter  component.  16 

7 Moderate  UHF  scintil  lation— focus  component.  17 

8 Severe  UHF  scinti 11 ation— probabi 1 i ty  densities  for 

real  and  imaginary  parts  of  scatter  component.  18 

9 Moderate  UHF  scinti  1 1 ation— probabil  i ty  densities  for 

real  and  imaginary  parts  of  scatter  component.  19 

10  Severe  UHF  sci nti 1 1 ation— probabi 1 i ty  densities  for 

phase  and  log-amplitude  of  the  focus  component.  20 

11  Moderate  UHF  scintillation— probabil ity  densities  for 

phase  and  log-amplitude  of  the  focus  component.  21 

12  Multiple  phase  screen  realization  of  the  severely 

scintillating  UHF  signal.  28 


3 


f 


FIGURE 

PAGE 

13 

Multiple  phase  screen  realization  of 
scintillating  UHF  signal. 

the  moderately 

29 

14 

Multiple  phase  screen  realization  of 
using  ionospheric  parameters  derived 
scintillating  UHF  signal. 

the  L-band  signal 
from  the  severely 

31 

15 

Simulated  severely  scintillating  UHF 
2-component  model . 

signal  using  the 

35 

16 

Simulated  moderately  scintillating  UHF  signal  using  the 
2-component  model . 

36 

r 


**r~~ 


'*■  •—  L£. 


SECTION  1 
INTRODUCTION 


Satellites  have  assumed  an  important  position  in  the  overall 
communication  capability  supporting  both  civil  and  military  operations. 
Consequently,  transionospheric  propagation  considerations  have  become  a 
necessary  part  of  system  design.  Experience  has  shown  that  the  ionosphere 
cannot  be  considered  a transparent,  homogeneous  propagation  medium,  even 
at  gigahertz  frequencies. 2 ‘‘Design  of  communication  signals  and  signal 
processors  requires  knowledge  of  the  signal  distortion  resulting  from 
transionospheric  propagation.  Models  of  the  ionospheric  distortion  process 
are  needed  for  these  designs  and  to  test  the  results  of  propagation  on 
specific  signals  and  systems.  A further  use  of  models  is  to  serve  as  a 
basis  for  extrapolation  from  natural  ionospheric  effects  to  effects  due 
to  ionospheric  perturbations  caused  by  high  altitude  nuclear  explosions. 
Since  little  data  on  nuclear  induced  scintillation  is  available,  propagation 
models  must  be  derived  on  a theoretical  basis.  These  models  should,  as  a 
minimum,  properly  predict  propagation  effects  through  a natural  ionosphere, 
given  an  adequate  description  of  the  ionosphere.  The  development  of  trans- 
ionospheric propagation  models,  then,  leads  to  derivation  of  models  for 
use  in  design  and  simulation  of  communication  systems  and  to  extraction  of 
pertinent  ionospheric  parameters. 

The  purpose  of  this  study  is  to  compare  the  results  of  two  models 
for  transionospheric  signal  scintillation.  The  basis  for  comparison  is  data 
taken  from  the  DNA  Wideband  satellite  experiment.1  The  models  used  are  the 
2-component  model5  developed  at  SRI  International  (SRI I),  and  the  multiple 
phase  screen  (MRS)  numerical  Fourier  propagation  simulation  model . 6 ’ 1 °For 


purposes  of  comparison,  two  time  periods  were  selected  from  Wideband  data 
taken  at  Ancon,  Peru,  on  December  16,  1976,  and  kindly  furnished  by 
Mr.  R.  C.  Livingston  of  SRII.  Selections  were  made  on  the  basis  of  the 
scintillation  index.  The  first  time  period  extends  over  a period  of 
40  seconds,  starting  at  4 hours,  47  minutes,  32  seconds  (UT).*  During 
this  period  signal  scintillation  was  severe.  The  scintillation  index  at 
413  MHz  was  1.00,  and  at  1239  MHz  the  scintillation  index  was  0.45.  During 
this  period  only  the  last  30  seconds  of  1239  MHz  data  was  used  because  the 
L-band  signal  was  used  for  phase  reference  during  the  first  10  seconds. 

The  second  time  period,  taken  from  the  same  satellite  pass, 
extends  from  4 hours,  49  minutes,  48.6  seconds,  to  4 hours,  50  minutes, 

28.6  seconds,  and  illustrates  a period  of  moderate  signal  scintillation. 
During  this  period  the  scintillation  index  for  the  UHF  signal  was  0.56. 


* Local  time  at  Ancon  is  five  hours  earlier  than  at  Greenwich.  Local 
time  was  23  hours,  47  minutes,  32  seconds,  December  15,  1976. 
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SECTION  2 

PARAMETER  ESTIMATION 


Experimental  data  is  recorded  as  quadrature  components  of  the 
received  signal.  The  Wideband  beacon  provides  an  S-band  signal  (2891.2  MHz) 
for  use  as  a phase  reference  for  coherent  detection  of  the  lower  frequency 
signals.  The  satellite  signal  received  on  the  ground  would  vary  in  ampli- 
tude and  phase  as  a result  of  changing  slant  range,  antenna  gain  pattern, 
and  total  electron  content  (TEC)  along  the  propagation  path,  even  if  the 
ionosphere  were  unperturbed.  As  the  satellite  progresses  through  the  pass, 
the  line  from  satellite  to  receiver  scans  a portion  of  the  ionosphere. 

The  process  of  scanning  translates  spatial  ionospheric  variations  into  tem- 
poral variations  of  the  received  signal.  The  variations  discussed  above, 
resulting  from  conditions  of  the  experiment,  take  place  slowly  and  are 
termed  trends.  The  initial  processing  step  is  to  remove  the  trends.  In 
Reference  5 , Fremouw  and  Rino  discuss  the  detrending  process  and  adopt  an 
upper  frequency  of  0.1  llz  for  trend  components.  This  limit  was  found  by 
use  of  the  criterion  that  detrending  should  not  reduce  the  S.  index  by  more 
than  "a  few"  percent.  It  is  clear  that  this  frequency  limit  is  dependent 
on  the  geometry  of  the  propagation  path.*  It  is  also  very  likely  that  it 
should  depend  on  ionospheric  structure.  For  this  work  we  have  used  the 
0.1  Hz  cutoff,  6 pole  Butterworth  numerical  filter  adopted  by  SRII.  Signal 
level  variations  with  periods  greater  than  ten  seconds  are  normally  of 
little  concern  in  operational  communication  systems. 


* Reference  5 reports  detrending  ATS-6  data  with  filter  time  constants 
sixty  times  as  long  as  those  used  for  Wideband  and  Transit  data. 


The  Wideband  satellite  data  is  presented  in  terms  of  quadrature 
components.  The  phase  of  the  signal  derived  from  the  quadrature  components 
is  made  continuous  by  limiting  the  phase  change  between  consecutive  samples 
to  plus  or  minus  tt  radians.  The  intensity  of  the  signal  is  the  sum  of  the 
squares  of  the  component  amplitudes.  Phase  and  logarithm  of  intensity  are 
low-pass  filtered  separately.  The  filter  output  is  then  subtracted  from 
the  input.  In  this  manner  log-intensity  and  phase  variations  with  periods 
greater  than  ten  seconds  are  removed  from  the  raw  data.  The  subtraction 
process  normalizes  the  intensity  by  the  average  over  the  filter  time 
constant . 


Detrended  amplitude  and  phase  for  the  Wideband  signals  at  415  MHz 
and  at  1239  MHz  are  shown  in  Figures  1 through  3.  Figure  1 shows  the 
detrended  intensity  and  phase  for  the  earlier,  severely  scintillating,  UHF 
signal  and  Figure  2 presents  the  moderately  scintillating  UHF  signal  during 
the  40  second  period  about  two  minutes  later.  Figure  3 shows  the  detrended 
L-band  signal  during  the  time  period  tvhen  the  UHF  signal  is  severely  disturbed. 
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Figure  1.  Severely  scintillating  UHF  signal —(detrended  intensity 
and  phase  from  Wideband  satellite  pass  at  Ancon,  Peru, 
on  16  December  1976). 
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Figure  2.  Moderately  scintillating  UHF  signal —(detrended  intensity 
and  phase  from  Wideband  satellite  pass  at  Ancon,  Peru, 
on  16  December  1976). 
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figure  3.  L-band  signal  during  severely  scintillating  UHF  conditions— 


(detrended  intensity  and  phase  from  Wideband  satellite  pass 
at  Ancon.  Peru,  on  16  December  1976). 


SECTION  3 

2-COMPONENT  MODEL 


The  2-component  signal  model  has  the  form 

E = EsEf  = (p+x+jy)exp(n+x+j4>)  , [1) 

where  p and  q are  real  constants,  x and  y are  joint  gaussian  zero 
mean  random  variables,  and  x anc^  ate  another  pair  of  joint  gaussian 
zero  mean  random  variables  statistically  independent  of  x and  y.  In 
this  2-component  model  scintillation  is  represented  by  two  statistically 
independent  multiplicative  random  components:  the  "scatter"  component  Eg, 
and  the  "focus"  component,  E^.. 

The  scatter  component  represents  the  rapidly  varying  signal  ampli- 
tude and  phase  caused  by  fairly  small-scale  ionization  structure.  It  is 
modeled  by  complex  gaussian  statistics  wherein  p is  the  unscattered  or 
specular  component,  and  x and  y are  a pair  of  correlated  zero-mean 
gaussian  random  variables,  x and  y need  not  have  equal  variances.  This 
scatter  component  is  associated  with  what  is  sometimes  referred  to  as  classical 
scintillation. 

The  focus  component,  E^.,  represents  the  more  slowly  varying 
components  of  amplitude  and  phase  caused  by  fairly  large-scale  ionization 
structure.  It  is  modeled  by  complex  log-normal  statistics  where  x and 
<J>  arc  correlated  zero-mean  random  variables  with  unequal  variances.  T’-e 
amplitude  fluctuations  associated  with  this  component  are  often  rather 
small;  the  primary  effect  of  this  component  is  usually  associated  with  the 
relatively  slowly  fluctuating  phase  of  the  composite  signal. 
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The  scatter  and  focus  components  of  the  detrended  signal  are 
separated  in  a second  filtering  process.  Because  this  model  assumes  statis- 
tical independence  of  the  two  components,  a 10-pole  Butterworth  numerical 
low-pass  filter  with  a cutoff  frequency  of  0.4  Hz  is  used  to  separate  the 
spectra  of  the  components.  The  output  of  the  low-pass  filter  is  taken  as 
the  focus  component.  Subtracting  the  focus  phase  and  log  intensity  from 
the  detrended  signal  gives  the  scatter  component.  Scatter  and  focus  compo- 
nents of  the  severely  and  moderately  scintillating  UHF  signals  are  given  in 
Figures  4 through  7. 

The  gaussian  nature  of  the  scatter  component  is  shown  in  Figures 
8 and  9.  In  these  figures  the  real  (p+x)  and  imaginary  (y)  parts  of 
the  complex  scatter  component  are  presented  as  a discrete  probability 
density  function.  The  normalization  performed  during  detrending  divides 
the  amplitude  terms  by  the  root  mean  square  of  the  total  receive^  signal 
amplitude.  A gaussian  probability  density  function  with  the  mean  and 
standard  deviation  calculated  from  the  data  for  each  case  is  plotted  with 
the  discrete  function  for  comparison.  We  note  that  the  imaginary  component 
has  an  essentially  zero  mean  for  both  time  segments  for  the  UHF  signal.  In 
the  severely  scintillating  UHF  signal  the  normalized  real  part  has  a mean 
of  0.28,  approaching  a Rayleigh  distribution.  The  mean  of  the  real  part 
of  the  scatter  component  for  the  moderately  scintillating  signal  is  nearly 
unity,  indicating  a largely  specular  process. 

Figures  10  and  11  show  the  distributions  of  the  amplitude  and 
phase  of  the  focus  component.  In  these  figures  the  normalized  value  of 
n + X is  the  logarithm  of  the  focus  component  amplitude.  The  phase  is 
the  quantity  (p  from  Equation  1.  The  discrete  distribution  functions 
plotted  in  these  figures  suffer  from  a limited  amount  of  data.  The  focus 
spectral  components  have  periods  ranging  from  2.5  to  10  seconds.  Therefore, 
the  40  seconds  of  data  used  here  contain  4 to  16  cycles  of  these  components. 
The  spikes  in  the  density  functions  clearly  match  the  maxima  and  minima 
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Figure  6.  Moderate  UHF  scintil  lation— scatter  component 
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Figure  8.  Severe  UHF  scinti 1 lation— probabil i ty  densities  for 
real  and  imaginary  parts  of  scatter  component. 
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Figure  9.  Moderate  UHF  scintillation— probability  densities  for 
real  and  imaginary  parts  of  scatter  component. 
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Figure  10.  Severe  UHF  scintillation— probability  densities  for 
phase  and  log-amplitude  of  the  focus  component. 
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Figure  1 


. Moderate  UHF  scintillation— probability  densities  for 
phase  and  log-amplitude  of  the  focus  component. 


of  the  corresponding  time  functions.  It  appears  that  the  discrete  density 
functions  would  as  likely  approach  gaussian  distributions  with  zero  means 
as  any  other  distribution  if  more  data  were  included  in  the  calculation. 
This  is  true  for  both  segments  of  the  UHF  signal. 


Statistics  for  the  detrended,  normalized,  continuous  phase  UHF 
signals  during  the  two  40  second  time  periods  are  given  in  Tables  1 and  2. 
Each  set  of  statistics  is  a result  of  processing  20,000  samples  of  the 
Wideband  satellite  signal.  The  parameters  in  the  table  are  defined  below 
for  reference.  Angle  brackets,  <>,  indicate  averages. 

- 2 2 1/2 

= Scintillation  index  = [(<I^>  - <I>  )/<I>  ] 
where  I is  intensity  = |e|“ 


a 

x 


= Root  mean  square  of 


= [<x  > 


2 

<x>  J 


1/2 


values  of 


x 


a = Root  mean  square  of  values  of  y . 

p = Normalized  coefficient  of  correlation  of  x and  y 
xy 

= <(x-<x>) (y-<y>)>/o  a r 

T = Time  interval  such  that  the  correlation  between 
x 

x(t)  and  x(t+Tj  = 1/e  = 0.37 
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Table  1.  Statistics  for  severely  scintillating  UHF  signal. 


Table  2.  Statistics  for  moderately  scintillating  UHF  signal. 
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SECTION  4 

MULTIPLE  PHASE  SCREEN  SIMULATION 

The  multiple  phase  screen  (MPS)  propagation  simulation  is  an 
analytical/numerical  technique  which  provides  a numerical  solution  for  the 
propagation  of  a plane  wave  through  a disturbed  ionosphere.  By  modeling 
the  ionosphere  as  a series  of  random  phase  screens  with  a power-law  power 
spectral  density  and  utilizing  the  actual  Wideband  satellite  geometry,  we 
attempt  to  "match1'  the  experimental  data  (in  a qualitative  sense)  and  thus 
obtain  some  information  on  possible  ionospheric  parameters. 

The  multiple  phase-screen  propagation  simulation  represents  the 
disturbed  region  by  a number  of  phase-screens  (10  for  this  work)  located 
in  the  disturbed  region  between  the  satellite  and  the  receiver.  Random 
phase  fluctuations  in  each  screen  are  generated  using  the  statistical 
properties  of  the  electron-density  fluctuations  as  determined  by  the 
electron-density  power  spectral  density.  A wave  (initially  plane  as  it 
enters  the  disturbed  region)  is  then  propagated  numerically  from  one  screen 
to  the  next  by  use  of  the  Fresnel-Kirchhoff  integral  equation  until  a solu- 
tion is  obtained  for  the  complex  electric  field  in  the  receiver  plane.  This 
technique  is  equivalent  to  a solution  of  the  parabolic  wave  equation  and  is 
thus  able  to  account  for  multiple  scattering.8  Since  the  phase-screens  are 
random,  the  signal  propagated  to  the  receiver  is  random  and,  if  desired, 
statistics  may  be  obtained  by  averaging  a number  of  different  simulations, 
each  based  on  a different  sequence  of  random  numbers. 

The  electron-density  power  spectral  density  is  a simulation 
input,  but  here  we  assume  a one-dimensional  power  spectral  density  of 
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the  form 


2 

a..  L 

scka)  = -N 

TT(1+KJL0“) 

where  is  the  standard  deviation  of  electron-density  fluctuation  and 

Lq  is  the  outer  scale  size.  Both  parameters  are  simulation  inputs.  Kj 
is  the  wavenumber  perpendicular  to  the  magnetic  field. 

We  assumed  a 100  km  thick  ionosphere  centered  at  an  altitude  of 
350  km.  To  model  the  actual  propagation  geometry  for  the  link  from  the 
Wideband  satellite  to  Ancon  we  used  constant  elevation  angles  of  15°  and 
30°  for  the  first  and  second  time  periods,  respectively.  The  actual 
satellite  elevation  angles  were  changing  slowly  during  the  40  second  time 
intervals  and  the  use  of  constant  angles  simplifies  the  analysis.  In  this 
manner  we  modeled  the  early  time  geometry  by  a layer  of  irregularities 
386  km  thick  with  the  center  1352  km  from  the  receiver.  In  the  late  time 
period  the  layer  of  irregularities  was  200  km  thick  centered  700  km  from 
the  receiver. 


The  best  match  of  MPS  calculated  data  to  the  actual  Wideband 
satellite  experimental  data  was  obtained  by  a visual  examination  of  a 
small  number  of  plots  of  intensity  and  phase  generated  by  the  MPS  propaga- 
tion simulation.  For  both  time  periods  we  found  that  an  outer  scale  size 

of  500  meters  gave  reasonable  fits  to  the  data.  During  the  first  time 

4 -5 

period  an  electron  density  standard  deviation  of  9.9  x 10  cm  was  used; 

during  the  second  time  period  an  electron  density  standard  deviation  of 
4-5  . 

5.3  x 10  cm  was  found  to  give  a good  tit  to  the  experimental  data. 

Values  for  both  the  outer  scale  size  and  the  electron  density  standard 
deviation  are  within  ranges  regularly  observed  in  the  equatorial  ionosphere. 

Since  the  propagation  simulation  results  arc  functions  of  distance 
in  a plane  normal  to  the  magnetic  field  direction,  the  simulation  results 
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may  be  converted  to  functions  of  time  by  dividing  the  simulation  computed 
distances  by  V . V is  defined  as  the  component  of  the  velocity  at  the 
ionospheric  penetration  point  of  the  propagation  path  perpendicular  to  the 
local  magnetic  field  direction. 

For  the  strong  U1IF  scintillation  case  a value  of  V of  780  m/sec 
was  required  to  match  the  observed  signal  intensity  decorrelation  time.  A 
value  of  V of  720  m/sec  was  required  to  match  the  intensity  decorrelation 
time  observed  for  the  moderate  UHF  scintillation.  The  geometry  of  the 
satellite  pass  indicates  that  a value  of  V of  about  500-700  m/sec  results 
from  satellite  motion  perpendicular  to  the  magnetic  field  direction.  Another 
100-200  m/sec  could  be  added  from  mean  atmospheric  drift.2  These  considera- 
tions indicate  a likely  range  of  V from  600-900  m/sec.  The  values  used 
here  are  within  this  range.  The  purpose  of  the  comparison  is  to  investigate 
adequacy  of  the  statistical  signal  description  in  the  two  models  — thus, 
correspondence  of  subsidary  quantities  such  as  Y , while  encouraging 
indicators  of  validity  for  MPS,  are  not  necessary  to  a valid  comparison. 

Note  that  the  simulation  parameters  which  match  the  Wideband 
experimental  data  are  not  unique.  They  were  chosen  simply  because  they 
yielded  the  closest  match  to  the  observed  signal  intensity  and  phase  from 
a small  set  of  simulations.  A somewhat  better  match  may  be  obtained  for 
a different  electron-density  standard  deviation  and  outer  scale  size. 

However,  the  signals  already  obtained  seem  to  match  the  data  quite  well 
and  we  do  not  expect  that  better  matches  will  result  in  much  change  in 
either  electron-density  standard  deviation  or  outer  scale  size. 

The  magnitude  and  phase  of  the  MPS  realizations  for  the  early 
time  period  are  shown  in  Figure  12.  These  can  be  compared  with  the  cor- 
responding detrended  Wideband  satellite  data  in  Figure  1.  Similarly,  the 
MPS  generated  amplitude  and  phase  for  the  late  time  period  are  presented 
in  Figure  15  for  comparison  with  Figure  2.  In  addition,  the  selected 
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Figure  12.  Multiple  phase  screen  realization  of  the  severely 
scintillating  UHF  signal.  (Compare  with  data  in 
Figure  1 . ) 
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Figure  13.  Multiple  phase  screen  realization  of  the  moderately 
scintillating  UHF  signal.  (Compare  with  data  in 
Figure  2.) 
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ionospheric  parameters  were  used  to  generate  a scintillating  signal  simula- 
tion at  L-band  for  the  early  time  period.  This  signal  is  shown  in  Figure  14 
for  comparison  with  Figure  3.  The  similarity  of  the  L-band  MPS  results  and 
the  detrended  Wideband  satellite  data  demonstrates  the  feasibility  of 
modeling  scintillation  effects  at  SHF  given  data  at  UHF.  L-band  propagation 
effects  for  the  late  time  period  were  not  modeled  because  the  relatively 
moderate  effects  at  UHF  indicate  that  effects  of  the  L-band  signal  pertur- 
bation will  be  small. 

Table  3 summarizes  statistical  parameters  indicating  severity  of 
scintillation  for  the  Wideband  data  and  for  signals  generated  using  the 
2-component  and  MPS  models. 
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Figure  14.  Multiple  phase  screen  realization  of  the  L-band  signal 
using  ionospheric  parameters  derived  from  the  severely 
scintillating  UHF  signal.  (Compare  with  data  in 
Figure  3.) 
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is  the  standard  deviation  of  phase  (radians). 


SECTION  5 

RECEIVER  SIMULATION 


To  investigate  the  effects  of  propagation  disturbance  modeling  on 

the  performance  of  system  models,  the  perturbed  signals  were  used  as  input 

to  a sampled  data  digital  simulation  of  the  DSCS-I1  AX/USC-28  modem.7  For 

this  purpose  a 75  bit  per  second  data  rate  was  used  with  coherent  detection 

of  a phase-shift  keyed  (PSK)  signal  with  differential  encoding  of  the  data. 

The  receiver  representation  used  a second  order  modified  Costas  phase-locked 

loop  with  a loop  bandwidth  of  21  Hz.  The  simulation  sampling  rate  was 

300/sec.  A mean  carrier-power-to-noise  density  of  35  db-Hz,  corresponding 

to  a mean  E,  /X  of  16  db  was  used  so  that  detection  errors  could  be  attri- 
b o 

buted  to  the  nature  of  the  signal  rather  than  to  receiver  noise.  Without 
signal  disturbances  the  error  rate  is  vanishingly  small  for  this  signal-to- 
noise  ratio.  The  receiver  model  can  accept  signal  input  data  in  each  of 
the  three  forms  discussed  here. 

The  form  of  the  direct  data  input  to  the  receiver  model  consists 
of  digitized  amplitude  and  phase  samples.  Since  the  Wideband  satellite  data 
is  sampled  at  a rate  of  500  samples/second  and  the  receiver  model  sample 
rate  is  500  samplcs/second,  an  interpolation  is  necessary.  Otherwise  the 
detrended  data  is  used  directly. 

Use  of  the  2-component  model  requires  tiiat  statistics  for  the 
scintillating  signal  be  derived.  Eight  parameters  are  required.  These 
are: 

• Scatter  scintil lat ion  index,  S ^ 

• Ratio  of  variances  of  \ and  v , a“/'~ 

x v 
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Correlation  coefficient  of  x and  y , p 

’ xy 

Decorrelation  time  (e  for  the  scattered  signal,  t_ 


Standard  deviation  of  the  log-normal  amplitude,  a 


X 


Standard  deviation  of  the  log-normal  phase,  0 


Correlation  coefficient  of  y and  $ , 


Xr 


Decorrelation  time  for  the  log-noi'mal  signal, 


These  quantities  are  given  in  Tables  1 and  2, 


The  random  scatter  signal  scintillation  components  are  computed 
at  each  sampling  time  using  a correlated  random  sampling  technique  applicable 
to  gaussian  random  processes.  The  correlated  sampling  technique  was  developed 
by  Hendrick  and  is  described  in  Reference  9.  A separate  set  of  correlated 
gaussian  random  variables  is  generated  for  use  in  the  representation  of  the 
focus  component  of  the  signal.  The  scintillating  signal  amplitudes  and 
phases  generated  by  this  process  are  presented  in  Figures  15  and  16.  The 
phase  as  plotted  in  Figure  15  is  constrained  to  an  interval  ranging  from 
-3n  to  3tt  radians..  Whenever  the  actual  phase  decreases  below  -3~  , 

2ir  is  added  to  all  subsequent  phase  values  and  plotting  continues.  The 
many  sharp  discontinuities  shown  in  Figure  15  are  results  of  the  addition 
of  these  2-  shifts  and  affect  only  the  plotted  phase,  not  the  receiver. 


The  form  of  the  '•IPS  inputs  to  the  receiver  model  is  the  same  as 
that  for  experimental  data.  These  are  the  signals  presented  in  Figures  12 
through  14 . 
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Figure  15.  Simulated  severely  scintillating  UHF  signal  using  the 
2-component  model.  (Compare  with  Figure  1.) 
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Figure  16.  Simulated  moderately  scintillating  UHF  signal  using 
the  2-component  model.  (Compare  with  Figure  2.) 


SECTION  6 

SIMULATION  RESULTS 


A phase  shift  keyed  bit  stream  which  carries  the  ASCII-coded 
message  "Simulation  of  coherent  PSK  satellite  communications  during  signal 
disturbances."  is  used  here  in  the  receiver  model.  The  biphase  modulated 
carrier  is  modified  by  the  signal  amplitude  and  phase  perturbations  from 
Wideband  data,  MPS  calculations,  or  the  2-component  model.  The  simulated 
receiver  attempts  to  track  the  carrier  in  a Costas  loop,  adjust  for  ampli- 
tude variations  with  an  AGC  circuit,  and  coherently  detect  the  data. 

The  following  is  an  example  of  the  message  output  from  the 
receiver  simulation  using  the  severely  scintillating  UHF  Wideband  data 
samples . 


'J  INUNCTION  OF  COHERENT  PSK 
Si'JJL*TION  OF; COHAPENT  PSK 
SIMULATION  OF  C'-HERENT#P,r.| . 
C IHUDATI0N  F COHERE  T#P.-L 
PIi-..jLH  i IClN#f'IF  COHEPr.HT  F"-l‘ 
SIMUNAT!  N OF  COHERENT  PSK 


SATELLITE  communications 
SATELL I TE  COtli  1U J I CATI 0 . S 
SATELLITE  ICOGMUN ICAUI  NS 

katull ite  c -i  s vjh i catioi ■ s 

SATELL I TE  CQMMUN I CAT I ON : 

rO 


DU 

NG 

sigfal 

DISTURBANCE  / 

DUB 

NG 

SIG  AL 

JSTUPBAI 

. 

■ 

. 

G IGNAL 

d rsTUPBPi 

g 

L 

N 

S I GNAL 

DISTURBAI- 

' . 

DUP 

NG 

SIGNAL 

D 3 TUP  HI 

. 

Outputs  using  signals  derived  from  the  models  for  this  time 
period  look  much  the  same  with  a different  array  of  errors.  Decoding 
errors  are  caused  by  phase  tracking  errors  and  phase  slipping  in  the  Costas 
loop,  as  well  as  by  amplitude  fades  below  the  threshold  for  reliable  data 
detection.  A summary  of  binary  error  statistics  from  the  simulations  is 
given  in  Table  4. 


Table  4.  Summary  of  error  statistics. 


An  explicit  bit  by  bit  demodulation  process  is  simulated  in  the 
receiver  model.  In  addition,  a conditional  binary  error  probability  is 


calculated.  The  conditional  error  probability  is  a function  of  the  phase 
tracking  error  and  the  signal  amplitude  at  each  sampling  time.  The  cal- 
culated bit  error  rates  given  in  Table  4 are  computed  by  averaging  the 
conditional  binary  error  probabilities  over  the  simulation  time  interval. 

These  are  the  average  error  rates  that  would  be  expected  for  very  long 

-4 

samples.  When  the  average  error  probability  is  of  the  order  of  10  , 

the  expected  number  of  bit  errors  in  a 40  second  period,  for  a bit  rate 
of  75/sec,  is  less  than  one.  However,  because  the  decorrelation  time  is 
long  compared  to  a bit  period,  the  errors  are  clumped  in  time.  Longer 
simulations  would  be  required  to  obtain  good  statistics  for  demodulated 
bit  errors  in  the  cases  involving  moderate  scintillation. 


SECTION  7 
CONCLUSIONS 


t 


Models  for  representation  of  radio  signals  that  have  propagated 
through  a disturbed  ionosphere  are  useful  for  generating  input  signals  for 
computer  simulation  of  communi cation  systems.  One  of  the  models  used  here, 
the  joint  gaussian  plus  log  normal  2-component  model,  represents  the  signal 
through  measured  statistical  parameters.  The  second  model,  the  multiple 
phase  screen  model,  represents  the  signal  by  generating  amplitude  and  phase 
samples  based  on  physical  propagation  effects.  In  the  comparison  of  these 
models  presented  here,  a digital  simulation  of  a BPSK  demodulator  was  found 
to  give  similar  demodulation  performance  for  the  two  signal  models.  The 
binary  error  rates  resulting  from  use  of  the  models  were  comparable  to 
results  obtained  from  use  of  the  data  taken  from  the  DMA  Wideband  satellite 
experiment . 


. 
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In  this  work  the  MPS  simulations  which  were  obtained  to  "match" 
the  Wideband  satellite  data  were  chosen  on  the  basis  of  scintillation  index 
and  phase  standard  deviation  from  a smaii  set  of  possible  cases.  Future 
work  of  this  type  should  be  directed  toward  matching  the  intensity  and 
phase  power  spectral  densities  of  the  received  signal  to  obtain  a match 
of  both  first  and  second  order  statistics.  In  addition  to  being  a further 
test  of  how  well  the  models  match  experimental  data,  these  quantities  impose 
requirements  on  automatic  gain  control  and  phase  tracking  loop  bandwidths. 

The  MPS  technique  has  advantages  for  communication  link  simula- 
tions because  the  parameters  required  are  related  to  ionospheric  structure 
rather  than  a priori  assumptions  of  signal  statistics. 


•ft-'” 
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ATTN:  CAPT  J.  Burton 
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DEPARTMENT  OF  THE  NAVY  (Continued) 
Fla va  1 Research  Laboratory 


DEPARTMENT  OF  THE  AIR  FORCE  (Continued) 
Deputy  Chief  of  Staff 


ATTN 

Code  7550 

J.  Davis 

Programs  & Analyses 

ATTN 

Code  5460 

Department 

of  the  Air  Force 

ATTN 

Code  7555 

ATTN 

PACSC , Maj  Paul 

ATTN 

Code  7580 

ATTN 

Code  6701 

J.  Brown 

Electronic 

Systems  Division, 

ATTN 

Code  6700 

T.  Coffey 

ATTN 

J.  Whelan 

ATTN 

Code  7500 

B.  Wald 

ATTN 

DCKC,  J.  Clark 

ATTN 

XRW,  J.  Deas 

Surface  Weapons 

Center 

ATTN 

YSEA 

ATTN 

Code  F31 

Naval  Surface  Weapons  Center 

Dahlgren  Laboratory 

ATTN:  Code  DF-14,  R.  Butler 

Navy  Space  Systems  Activity 
ATTN:  A.  Hazzard 
ATTN:  Code  52 

Office  of  Naval  Research 
ATTN:  Code  461 
ATTN:  Code  421 
ATTN:  Code  402 
ATTN:  Code  420 

Office  of  the  Chief  of  Naval  Operations 
ATTN:  0p-604 
ATTN:  Op-941 

Strategic  Systems  Project  Office 

Department  of  the  Navy 

ATTN:  NSSP-2722,  F.  Wimberly 
ATTN:  NSP-2141 

DEPARTMENT  OF  THE  AIR  FORCE 

Aerospace  Defense  Command/DC 

Department  of  the  Air  Force 
ATTN:  DC,  Mr.  Long 

Aerospace  Defense  Command/XPD 

Department  of  the  Air  Force 
ATTN:  XPOOQ 
ATTN:  XP 

Air  Force  Avionics  Laboratory 
ATTN:  AAD,  W.  Hunt 
ATTN:  AAD,  A.  Johnson 

Air  Force  Geophysics  laboratory 
ATTN:  OPR-1 , J.  Ulwick 
ATTN:  PHP,  j.  Aarons 
ATTN:  LKB,  K.  Champion 
ATTN:  OPR,  A.  Stair 
ATTN:  PHD,  J.  P.uchau 
ATTN:  PHD,  J.  Mullen 

Air  Force  Technical  Applications  tenter 
ATTN:  TF/Maj  Wiley 

Air  Force  Weapons  Laboratory 
ATTN:  SHI 
ATTN:  DYC,  .1.  Kanin 
ATTN:  DYF , J.  Earry 

Deputy  Chief  of  Staff 

Research,  Development,  A Acq. 

Department  of  the  Air  Force 
ATTN:  AFRDO 


Foreign  Technology  Division,  AFSC 
ATTN:  FTDP,  B.  Ballard 
ATTN:  N1CD,  Library 

Rome  Air  Development  Center,  AFSC 
ATTN:  Documents  Library/TSLD 
ATTN:  OCSE,  V.  Coyne 

Rome  Air  Development  Center,  AFSC 
ATTN:  EEP 

Space  and  Missile  Systems  Organization/SK 
Air  Force  Systems  Command 
ATTN:  SKA,  M.  Clavin 

Strategic  Air  Conmand/XPFS 
Department  of  the  Air  Force 

ATTN:  DOK,  Chief  Scientist 
ATTN:  NR 7 

DEPA  TMJ  NT  Nl  ■ 

Department  of  Energy 
Alhuguergup  Operations  Office 

ATTN:  Dui.  Con.  for  D.  Sherwood 

Department  of  Energy 
library  Room  G-042 

A 1 TN : oc . I . • ( I 

Office  of  Military  Appl’intion 
Department  of  Energy 

ATTN:  Poc.  Con.  for  1 . dale 

iOVEI  NM  N 

Central  Intelligence  Agency 

ATTN:  RD”  I , Rm.  * G4!  . Hq.  nidg. 
for  OS] f O'  TO 

Department  of  Commerce 
National  P.ureau  of  Standard' 

ATTN:  R,  Moore 

Department  of  Transportation 
Office  of  the  Secretarv 
ATTN:  R.  Doherty 
ATTN:  R.  Lewis 

Institute  for  Telecommunications  Sciences 
National  leleconmuni cation'  1 Info.  Admin. 
ATTN:  W.  Iltlaut 
ATTN:  0.  frontiie 
ATTN:  L.  Perry 
ATTN:  A.  Dean 
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''THER  GOVERNMENT  AGENCIES  (Continued) 


DEPARTMENT  0£  DEFENSE  CONTRACTORS  XConti nued ) 


National  Oceanic  & Atmospheric  Admin. 
Environmental  Research  Laboratory 


General  Electric  Company 
ATTN:  E.  Reibert 


Department 

of  Commerce 

ATTN 

G.  Reid 

General  Electric  Co. -TEMPO 

ATTN 

R.  Grubb 

Center  for 

Advanced  Studies 

ATTN 

W.  Knapp 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

ATTN 

DAS 1 AC 

ATTN 

T.  Stevens 

Aerospace  Corporation 

ATTN 

D.  Chandler 

ATTN 

I.  Garfunkel 

ATTN 

M.  Stanton 

ATTN 

SMEA  for  PWW 

ATTN 

S.  Bower 

General  Electric  Tech.  Services  Co.,  Inc. 

ATTN 

F.  Morse 

ATTN 

G.  Mi  liman 

ATTN 

V.  Josephson 

ATTN 

T.  Salmi 

General  Research  Corporation 

ATTN 

J.  Carter 

Santa  Barbara  Division 

ATTN 

N.  Stockwell 

ATTN 

J.  Ise,  Jr. 

ATTN 

D.  Olsen 

ATTN 

J.  Garbarino 

Analytical 

Systems  Engineering  Corp. 

Geophysica 

Insti tute 

ATTN 

Radio  Sciences 

University 

of  Alaska 

ATTN 

N.  Brown 

Berkeley  Research  Associates,  Inc. 

ATTN 

T.  Davis 

ATTN 

J.  Workman 

ATTN 

Technical  Library 

Boeing  Company 

GTE  Syl vania,  Inc. 

ATTN 

G.  Mall 

Electronics  Systems  Grp. -Eastern  Div. 

ATTN 

D.  Murray 

ATTN 

M.  Cross 

ATTN 

G.  Keister 

ATTN 

J.  Kenney 

HSS,  Inc. 

ATTN 

D.  Hansen 

Brown  Engineering  Company,  Inc. 

ATTN 

R.  Deliberis 

University 

of  Illinois 

Department 

of  Electrical  Enaineering 

Universi ty 

of  California  at  San  Diego 

ATTN 

■ . Yeh 

ATTN 

H.  Booker 

Insti tute 

for  Defense  Analyses 

Charles  Stark  Draper  Lab.,  Inc. 

ATTN 

. Aein 

ATTN 

J.  Gilmore 

ATTN 

'.  Benqston 

ATTN 

D.  Cox 

ATT!. 

f.  Bauer 

A ’ T N 

' . .-.'olfhard 

Computer  Sciences  Corporation 

ATTN 

H.  Blank 

Interna*  ;o*'al  Tel.  \ Telegraph  "'orporation 

ATTN 

Technical  Library 

COMSAT  Labs. 

ATTN 

G.  Hyde 

JAYCOR 

ATTN 

R.  Taur 

ATTN 

S.  ioldran 

Cornell  University 

Johns  • r:  ► v niversitv 

Department  of  Electrical  Engineerinn 

Applied  Fh 

rics  tab. 

ATTN 

D.  Farley,  Jr. 

A TTN 

Document  Librarian 

ATT*, 

’ . Tote- ra 

Electrospace  Systems,  Inc. 

ATTN 

H.  Logston 

Kaman  Sciences  rcrpp»*a t i e* 

ATTN 

T.  teacher 

ESL,  Inc. 

ATTN 

J.  Marshall 

Lawrence  Livermore  Laborat  * . 

ATTN 

C.  Prettie 

I Ini  versi  ty 

of  Cali fomia 

ATTN 

J.  Roberts 

AT  TN 

c.  • . • ■ • . ' ■ . 

ATTN 

• . • • -46,  . • • 

Ford  Aerospace  A Communl cat  ions  Corporation 

ATTN 

J.  Mattinqley 

Linkabi t Corporati  'r 

AT  T*. 

’.  acobs 

General  Electric  Company 

Space  Division 

LocFheed  Vissile-  *»  Sr-ao-  . 

ATTN 

M.  Eortner 

ATTN 

Dept.  60-12 

ATTN:  T.  Churchill 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued; 


Lockheed  Missiles  and  Space  Co.,  Inc. 
ATTN:  W.  Imhof 
ATTN:  M.  Walt 
ATTN:  R.  Johnson 


Sandia  Laboratories 
1 ivermore  Laboratory 

ATTN:  Doc.  Con.  for  B.  Murphey 
ATTN:  Doc.  Con.  for  T.  Cook 


Los  Alamos  Scientific  Laboratory 


Sandia  Laboratories 


ATTN: 

Doc. 

Con 

ATTN: 

Doc. 

Con 

ATTN: 

Doc. 

Con 

M.I.T.  Lincoln  Laboratory 


McDonnell  Douglas  Corporation 


ATTN: 

Doc. 

Con. 

for 

Space  Project 

ATTN: 

Doc. 

Con. 

for 

D.  Thornbrouah 

ATTN: 

Doc. 

Con. 

for 

D.  Dahlqren 

ATTN: 

Doc. 

Con. 

for 

3141 

ATTN: 

Doc. 

Con, 

for 

W.  Brown 

ATTN: 

D.  Towle 

ATTN: 

L.  Lough! in 

Science  Applications,  Inc 

ATTN: 

F.  Waldron 

ATTN:  D.  Sachs 

ATTN: 

C.  Clark 

ATTN:  L.  Linson 

ATTN:  D.  Hamlin 


ATTN 

W.  Olson 

Science  Applications,  Inc. 

ATTN 

G . Mroz 

Huntsville 

Division 

ATTN 

J.  Moule 

ATTN 

D.  Divis 

ATTN 

N.  Harris 

Science  Applications,  Inc. 

Mission  Research  Corporation 

ATTf! 

SZ 

ATTN 

R.  Scott 

ATTN 

D.  Knepp 

SPI  International 

ATTN 

M.  Scheibe 

ATTf, 

G.  Price 

ATTN 

D.  Sowle 

ATTN 

A.  Burns 

ATTN 

F.  fajen 

ATTN 

G.  Smith 

ATTN 

R.  Bogusch 

ATTN 

C.  Pino 

ATTN 

R.  Hendrick 

ATTN 

D.  Neilson 

ATTN 

S.  Gutsche 

ATTN 

W.  Chesnut 

5 cy  ATTN 

Document  Control 

ATTN 

G.  Carpenter 

ATT*' 

Faron 

Mitre  Corporation 

ft  T T • ■ 

V.N  Jaye 

ATTN 

C.  Callahan 

ATTN 

P.  Livingston 

ATTN 

G.  Harding 

p.  Leadabrand 

Mitre  Corporation 

Tri-Con',  Inc. 

ATTN 

M.  Horrocks 

f.  7 T !. 

: D.  Murray 

ATTN 

VI.  Foster 

ATTN 

W.  Hall 

U i efens(  S ! pa<  e Sys.  in  uj 

ATT!; 

S • Al  t huler 

Pacific-Sierra  Research  Corporation 

T T*i 

ATTN 

L.  Field,  Jr. 

ATT*; 

: r . Dee 

Pennsylvania  State  University 

"tah  ftate 

! 'ni  v»".  l ty 

ATTN 

Ionospheric  Research  Lah. 

Center  for 

' esearc'-  • » ' • ' 

Space  Science  I ahr  rat ■ - 

Photometries,  Inc. 

f.  7T’, 

: '• . Bai  er 

ATTN 

I.  Kofsky 

: L.  Jensen 

Physical  Dynamics,  Inc. 

Vf Si  : ■ ■ • . 

Inc. 

ATTN 

E . Fremouw 

rarpenter 

? t D Associates 

Riverside 

Research  Institute 

ATTf 

C.  MacDonald 

AT'*. 

: V,  Trapani 

ATT’ 

F.  Gilmore 

ATTf 

W.  Wright,  Jr. 

ATTf 

W,  Yarza s 

ATTf 

P.  Lelevier 

ATTf 

R.  Turco 

ATTf 

H.  Dry 

ATTf 

B.  Gabbard 

Rand  Corporation 

ATTN 

: C,  Crain 

ATTN 

: T,  P.edrozian 
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